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disturbance  signals.  In  tl\c  SOAS , p-  reappears  as  a significant 

factor  in  the  quasilincar  requirements  which  may,  depending  on 

the  numbers  involved,  completely  vitiate  its  banishment  as  an 

uncertainty  factor.  The  development  of  a quantitative  design 

theory  for  t)io  SOAS  pinpoints  tlie  two-loop  SOAL  extension  needed 

to  overcome  this  SOAS  f.imitation,  and  the  development  of  a 

similar  SOAL  quantitative  design  theory.  In  the  latter  o 

the  ’ 

disappears  from  both^adapti ve  and  quasilinear  conditions, 

but  reappears  as  a factor  in  the  rate  of  adaptation  of  the 

^system.  It  may  be  banished  from  here  too,  by  means  of  a 

three-loop  self-oscillating  system  (SOANL),  for  which 

the  SOAIj  design  theory  is  applicable  with  minor  extensions. 
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A SYNTHESIS  THEORY  FOR  MULTIPLE-LOOP  OSCILLATING  ADAPTIVE  SYSTEMS 

Isaac  Horowitz'  and  Aharon  Shapiro 


ABSTRACT 

The  multiple-loop  self-oscillating  adaptive  (SOAL)  system 
is  presented  as  a natural,  logical  means  of  overcoming  a serious 
limitation  of  the  single-loop  self-oscillating  system  (SOAS) . 

Both  structures  have  the  property  of  zero  sensitivity  to  plant 
high-frequency  gain  uncertainty  p = ^max'^^min  ’ factor 

which  is  generally  responsible  for  large  'cost  of  feedback’.  It 
is  however  necessary  to  design  these  systems  such  that  the  response 
is  essentially  quasilinear  to  the  desired  clar^  of  command  and 
disturbance  signals.  In  the  SOAS,  p rea^  as  a significant 

factor  in  the  quasilinear  requirements  whicii  may,  depending  on 
the  numbers  involved,  completely  vitiate  its  banishment  as  an 
uncertainty  factor.  The  development  of  a quantitative  design 
theory  for  the  SOAS  pinpoints  the  two-loop  SOAL  extension  needed 
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to  overcome  this  SOAS  limitation,  and  the  development  of  a 

similar  SOAL  quantitative  design  theory.  In  the  latter,  p 

the 

disappears  from  both^adaptive  and  quasilinear  conditions, 
but  reappears  as  a factor  in  the  rate  of  adaptation  of  the 
system.  It  may  be  banished  from  here  too,  by  means  of  a 
three-loop  self-oscillating  system  (SOANL)  , for  which 
the  SOAL  design  theory  is  applicable  with  minor  extensions. 


A SYNTHESIS  THEORY  FOR  MULTIPLE-LOOP  OSCILLATING 


ADAPTIVE  SYSTEMS 


I . INTRODUCTION 

It  is  well  known  that  high  frequency  dither  can  linear- 
ize the  response  of  a certain  class  of  nonlinear  elements  to 
slower  and  smaller  amplitude  signals  [1-3].  This  phenome- 
non is  explicable  by  multiple-input  describing  function 
theory  [4-7],  and  has  recently  received  more  abstract  and 
rigorous  treatment  [8].  The  adaptive  property  of  such  dith- 
er was  later  recognized  [9,  10  and  especially  11],  Struct- 
ures similar  to  dithered  systems  were  among  those  promoted 
in  the  USAF  sponsored  competition  for  adaptive  flight  con- 
trol systems  [12,  13].  Some  versions  were  among  the  suc- 
cessful ones,  used  for  example  in  the  X-15,  X-20  involving 
fantastic  parameter  uncertainty  factors  of  more  than  one 
thousand  [14-16],  These  systems  contain  strongly  nonlinear 
elements  but  their  system  properties  have  significant  linear 

aspects.  Thus,  aside  from  their  inherent  value,  they  are 

the  linear  and  nonlinear 

of  interest  as  a transition  between^adaptation  philoscphics . 

A variety  of  dithered  feedback  systems  exist.  In  most, 
the  amplitude  of  the  dither  is  used  as  a parameter  identi- 
fier, but  in  some  the  phase  is  also  used.  Single  and  tv/o- 


r 
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loop,  self-oscillating  and  externally  excited  systems  have 
been  described  [17-22].  In  all,  both  the  quasilinear  and 
the  adaptive  aspects  are  essential  features  which  must  be 
properly  designed  into  the  system.  The  quasilinear  aspect 
has  received  thorough  quantitative  treatment  [24],  but  the 
adaptive  aspect  to  a much  lesser  extent.  Since  parameter 
uncertainty  is  the  motivation  for  these  systems,  surely  the 
actual  extent  of  uncertainty  and  the  desired  system  response 
tolerances,  ought  to  appear  as  significant  design  parame- 
ters. It  is  reasonable  to  expect  that  an  optimum  design  to 
handle  an  uncertainty  factor  of  say  1000,  should  be  signi- 
ficantly different  from  one  for  an  uncertainty  factor  of 
100.  But  this  is  not  so  in  most  of  the  literature.  There 
have  so  far  appeared  quantitative  design  theories  only  for 
the  single-loop  self-oscillating  (SOAS)  and  the  single-loop 
externally  excited  (EEAS)  adaptive  systems  [25,  26].  This 
paper  presents  such  a design  theory  for  two  versions  of  the 
two-loop  self-oscillating  system  (SOAL  and  SOANE).  The 
SOAL  structure  used  here  is  not  new,  although  it  was  inde- 
pendently derived  in  a straightforward  logical  manner  (as 
will  be  seen) , in  the  process  of  trying  to  overcome  SOAS 
limitations.  The  quantitative  design  theory  is  new  and 
reveals  the  the  superiority  of  the  SOAL  structure  over  oth- 
er two-loop  versions.  The  quasilinear  feastures  common  to 
all  dithered  systems  are  first  reviewed. 

Review  of  Quasilinear  Theory  [ 24 ] 

Let  the  input  to  the  nonlinearity  N (Fig.  1 ) be 
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X = A sin  + x^(t)  - + x^(t)  (la,b) 

Let  0)^  be  the  bandwidth  of  X^Cjoi)  = s = jo)  and 

a,  P specific  suitable  numbers  (a  = 3 = 3 will  be  used 
in  subsequent  numerical  examples).  For  a certain  class  of 
nonlinearities j if 

(0 

(2a, b) 


A 


sup|x^(t) 1 ^ a > “f  ^ g 


then  y,  the  output  of  N may  be  written  with  good  accuracy 
(depending  on  the  values  of  a,  3 [24]) 

y = N^A  sin  m^t  + N^x^(t);  = ^ , ^ (3a-c) 

with  M a real  parameter  of  N.  The  ideal  relay  with  +B 

saturating  levels  will  be  always  used  henceforth  for  illus- 

trative  purposes,  for  which  M = — . Thus,  N presents 

the  gain  to  the  faster  and  larger  'oscillating  component 

X and  N-  = N /2  to  the  slower  and  smaller  forced  com- 
o f o 

ponent  x^  of  (la,b).  Accordingly  in  Fig.  1 the  forced 
component  system  response  functions  to  command  and  disturb- 


ance inputs  are 

A Lc  A 

Tf(s)  = -Y  " '^r+Tj  » "^d 


D 


1 + ’ 


(4a,b,c) 


^ Z 


dz 


d . -h 


1 + L. 


t>  -».  • i-f 

with  forced  and  oscillating  component  loop  transfer  func-.- 
tions 


1.^  = G-,G.,N^.P  = G,G^N^KP 


’1"2‘  f 


1"2  f ‘h 


L = G,G„N  P = ?L 
o 1 2 o 


N 


N. 


o _ 2B 
2 TiA 


f ’ 

(5a,b,c) 


-4- 


In  the  above  P = K is  the  constrained  plant  and  it 

is  assumed  in  the  meantime  that  the  only  uncertainty  is  in 
the  hif,h  frequency  plant  gain  factor  K , defined  by 


lim 

S->-oo 


PCs)  “ _1 
oP 


K 


(6) 


where 


is  the  excess  of  plant  poles  ovei"*  zeros. 


II.  REVIEW  OF  SO AS  ADAPTIVE  THEORY 

There  is  a logical  and  natural  transition  from  the 
SO  AS  to  the  two-loop  SOAL,  which  rendex^s  the  latter  theory 
quite  transparent  once  the  former  is  well  understood.  First, 
note  th£it  if  the  excess  of  poles  over  zeros  of  I.^(s)  of 

(5b)  is  5 3 , then  from  describing  function  theoi’y  ^nd  the 

infinite  gain  of  the  ideal  relay  (or  lai''ge  enough  gain  for 
other  nonlinear  element) , there  is  a limit  cycle  at  , 

defined  by 

Arg  Lo(3Wo)  = with  |L^(ja)^)I  = ^ 

(7a,  b) 

Since  all  elements  of  (7b)  except  K,  A are  fixed,  K/A 

must  be  a constant,  giving  zero  sensitivity  of  L^Cjoj)  = 

to  K uncertainty.  Incidentally,  the  rigorous  requiivments 
for  excellent  accuracy  of  the  describing,  function  [27]  arc 
easily  achieved  in  SOAS  and  SOAL  design.  Second,  in  £iny 
practical  system  there  must  be  a bound  m on  ttie  tolerable 
amplitude  of  the  limit  cycle  at  the  system  output,  i.e.  in 
Fig.  1 
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iGo  K (jw  )1 


= ^ Is  S 


I ^ m (8) 


Third,  the  quasilinearity  conditions  (2a,b)  compel 
specification  of  the  extreme  command  and  disturbance  signals 
for  which  quasilinearity  is  to  be  valid.  It  is  found  [25] 
that  these  are  the  inputs  giving  maximum  forced  component 
plant  output  at  in  Fig.  1,  and  denoted  by 

= Rg  Tjp  of  (4)  for  the  extreme  command  input,  and 
Zgj  = -D^  of  (4)  for  the  extreme  disturbance  input. 

The  corresponding  extreme  input  to  N is  denoted  by 

with  Z^^  (s)  = d(t)  , 

r(t)  never  or  rarely  occur  simultaneously  then  the  larger 
t each  w must  be  used.  If  the  two  extremes  do  occur 
together,  then  their  sum  must  be  used.  Fourth,  a good 
optimization  criterion  is  the  minimization  of  of  (7), 

because  the  bandwidth  of  is  thereby  minimized  [28]. 

In  the  SOAS,  conditions  (2a, 8)  usually  determine 

9 R 

^“o^min'  ’ follows.  Write  Z^Cjo))  = G2KPj^Xg(jm)  = 

2B 

— T Go  K . P,  X (jm)  because  A/K  is  invariant,  as 

ttA  . 2 min  h e 


previously  noted  following  (7).  Use  the  above  to  eliminate 
BG„P,  in  (8),  giving 


<3“o)  ^ HTT^  |2e'j“o>l= 
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The  final  important  constraint  is  (2a)  in  which  the  worst 

case  is  when  x and  A = A . . It  is  shown  in  [25’J 

f e min 

that 

that  the  best^can  be  done  results  in  (2a)  being  equivdilent 
to 


X 

e 

A . 

min 


( jm) 


1 

a oj  ’ 


o 


(10) 


The  importance  of  conditions  (9,  10)  is  illustrated  by  the 
following  desing  example. 

K 

Design  Example.  P = — , Kg  [1,100]  . The  extreme  command 
7 5 

input  R (s)  = — and  the  disturbances  are  thought  to  be 
o s 

negligibly  small.  The  maximum  tolerable  oscillation 
amplitude  m = 1 . The  desired  system  response  function  is 

T^(s)  = P;  Qi  ^ where  4>(s)  contains  the  far-off 

s'^  + O.  ls+0. 01 

poles  and  zeros  which  may  be  added,  if  worthwhile. 

Design.  In  the  above  Yi  = 2 K /m  K . = 200, 

s — '1  max  min  ’ 

Z = Z = R T^  and  the  resulting  YilZ  (jw)|  is 
s)<etched  in  Fig.  2 for  $(s)  = 1 as  is  1/a  w (a  = 3) 

of  (10),  as  if  testing  any  w if  it  may  be  used  for  . 

Conditions  (9,10)  dictate  that  (m^)  is  at  the  inter- 

action  of  the  two  curves,  with  value  20  rps  . Suppose 

$(s)  = — - — y is  acceptable.  The  new  a,|Z  (juOl  gives 
(s  + l)'^  ^ ^ 

u>  = 4.5  . Clearly,  it  is  desirable  to  add  as  many 

o • 

min 

far-off  poles  as  tolerable.  Condition  (2b)  is  then  satisfied 
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by  a large  margin. 

It  is  seen  that  y- 


in  determining  (m  ) . 

^ o min 


2K 


max 


mK 


is  decisively  important 


max 


It  is  interesting  that  the 
K 

max 


uncertainty  in  K , expressed  by  > which  v;as  banished 

min 

from  Lj  , returns  as  a constraint  which  is  very  important 
in  determining  ^“o^min  with  it,  the  bandvzidth  of  . 

If  is  large  enough,  the  resulting  bandwidth 

may  be  larger  than  that  required  by  a linear  time-invariant 
feedback  design  satisfying  the  same  quantitative  specifica- 
tions [28]. 


Limitations  due  to  Disturbances 


Disturbances  impose  a serious  limitation  on  the  SOAS. 
This  is  seen  by  asking:  Given  determined  by  as 

in  the  above  example,  for  what  class  of  disturbances  D 
in  Fig.  1,  is  quasilinearity  maintained?  Consider  then 


Z^Cjm  ) = Z,  (jm  ) = 
e o de  o 


-h  . 

T+L7  ^e  (jeo^)  = D^(ja)^),  because  from 


(5b,  7b)  = -0.5  . In  the  limiting  case  (from  9,  10) 

Y-,  |Z  , (jm  )|  = Y-,  |D  (jo)  )|=l/3cj  , which  determines  the 

tolerable  D class.  Thus  in  the  above  example  if 

D = k/s  , then  set  200  k /oj  = 1/3  m , giving 

e max  o o ' 

k = 1/600  - a very  small  tolerable  step  disturbance  - 

the  assumed  negligible  disturbance  may  not  be  ignorable 

k 


after  all.  If  D, 


(s  + a) 


Tj-  with  w < (1)  , then 

z o 


I 

I 


k . 1 w 

200  = 3—  gives  . Any  larger  disturbances 

o ° 

violate  quasilinearity  condition  and  the  response  to  them 
becomes  very  sensitive  to  the  relative  phases  of  the 
sinusoil  and  the  disturbance  input. 

In  the  above,  the  problem  of  design  for  specific 

disturbance  attenuation  was  not  even  considered.  Suppose 

there  are  significant  disturbances  w'^ich  must  be  attenuated. 

Whatever  disturbance  attenuation  design  philosophy  [28]  is 

used,  it  results  in  some  extreme  plant  output  for  the 

extreme  disturbance  input  . As  in  the  aibove , it  is 

necessary  that  ^ ‘ ~ ’ 

there  is  no  solution  if  y^^^  ^ 1/°^  • How'-  er  if  the  excess 

of  poles  over  zeros  of  D^(s)  > 2 , then  there  is  always  a 

solution,  but  (m  ) . may  then  be  exceedingly  large.  Thus, 

’ o min  o j 0 j 

2 . 2 

in  the  above  example,  = J/s  gives  ’ 

so  (w  ) . = 6000  J . A linear  time-invariant  design 

o min  ^ 

might  be  far  better  in  terms  of  loop-bandwith. 

Completion  of  Design  Example 

If  disturbance  attenuation  is  of  no  concern,  then  the 
design  is  basically  complete,  once  has  been  chosen  [25]. 

For  0)  < (1)  , ll.  (ja))|  should  be  chosen  >1  of  course, 

but  not  necessarily  very  large,  because  as  such  is  not 

needed  for  disturbance  attenuation  or  sensitivity  reduction. 

It  is  desirable  to  have  large  |d  Arg  I.^(iu)/dml  at  , 

so  that  the  value  of  ui  is  insensitive  to  changes  in  plant 
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dynamics.  Such  an  Arg  L^(jw)  is  associated  vyith  fast 


decrease  of  |L„(ja))|  near  w 
' O ‘ o 


, which  in  itself  is  also 
desirable.  If  disturbance  attenuation  is  of  concern  then 


Lj  = 0.5  must  be  shaped  to  obtain  the  desired 

disturbance  attenuation,  vjhich  must,  as  previously  shown. 


be  compatible  with  the  chosen  u 


In  fact,  from  the 


previous  discussion,  m due  to  quasilinearity 

min 


constraints  is  very  likely  much  larger  than  that  needed  for 
disturbance  attenuation  alone.  As  an  example  of  design  for 
(compatible)  disturbance  attenuation  suppose  one  wants  it 
(i.e,  |T^(ja))|  of  4b)  to  be  ^ B^(w)  , a straight  line 

on  a Bode  plot  (db  vs  log  w)  , with  B^(l)  =-5db  , 

B^(*l)  = -30db  . It  is  given  that  this  bound  applies  for 
01  < 1.2  . For  all  other  o)  it  is  required  that 


1+L. 


^ 2.3db  . See  [28,  29  Sec.  10  . 15]  how  such  bounds 

may  arise.  In  any  case,  no  matter  what  design  approach  is 
used,  disturbance  attenuation  is  a matter  of  small  enough 
|T^(j(i))|  vs  01  over  a sufficiently  large  oi  range. 
Different  techniques  only  give  different  B^(io)  for  T^Cjoj) 

on 


Since  - , the  B^(a))  generate  bounds  BCw) 


L^( jw)  , shown  in  Fig.  3,  which  together  with  the  previous 


constraint  of  a)_  = 4.5  and  L^Cju^)  = -0.5  (point 


o min 


A in  Fig.  3),  complete  the  specifications  on  l.jj,(j(j)  . 
In  this  case,  L^(ja))  is  easily  shape'  to  satisfy  the 


B(a3)  and  ai^  = 4 . 5 , as  shown  in  Fig.  3,  with 


(1.28)10°(s‘+30s+225) 


Lj-(s)  = 


s(s^  + 6.44s  + 21.2)  (r.^  + 15s  + 225)^ 


Ik 
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The  prefilter  is  obtained  from 


T^(s) 


[L^/(l+L^)]  (1+lOs+lOOs^) (1+s)^ 


Constraint  (10)  determines  G„(s)  by  use  of 


% " ^ giving  1 862(30))  1 ^ 


2 

TT  acj 


2K  . 
min 


Z^(ja)) 


NfG2K 


Pj^(  jm) 


for 


w 0)  , 

o ’ 


with  I I the  larger  of  ■>  | | , assuming 

simultaneous  occurrence  of  extreme  r(t)  and  d(t)  is 
very  rare.  In  this  case,  it  is  easily  found  that 


Z (jco) 
er 


for  u $ 0) 


o 


so 


= R 


IS 


used.  A satisfactory  choice  is 


BG, 


•76tt(1+10s)(s+3.22s  + 21.2) 


(1+lOs+lOOs^) (1+s)^ 


Gj^  is  found  from 


L - ^ 
^f  " ttA 


B can  be  chosen  for  convenif 
in  hardware  implement  at  ion . The  value  8=  ISir  was  here  chosen. 


A . = 0 . 5 is  knovm  from  knowledge  of 

min  1 


BGz  , Gj^ , 


KPv 


Motivation  for  the  Two-Loop  SOAL 


It  is  seen  from  the  above  that  the  presence  of 


K /K  . in  Yt  of  (9)  is  a very  serious  impediment  in 
max  min  1 


SOAS  design,  leading  to  ^‘^o^min  unnecessarily  large  for 


the  adaptive  and  disturbance  attenuation  needs  of  the 


system.  Suppose  K /K  . could  be  eliminated  with 
^ ^ max  rnin 


Yz  = Z/m  used  instead  of  Yq^  in  (9).  Then,  in  the  above 
example  in  Fig.  2,  the  intersection  of  Y2 1 ^ with 


l/3u)  gives  (o)  ) . = 3.3  , in  p3acc  of  4.5rps,  an 

o min  ' 


J 
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improvement  of  1.79  octaveG.  Note  thfit  constraint  (lb)  is 

easily  satisfied  even  at  this  value  of  u . But  in 

o 

general  there  is  some  possibility  that  (lb)  becomes 

significant,  in  the  SOAL  design.  For  D = k/s  or 

9 ^ K 

k/(s+a)  etc.,  tolerable  k is  exactly  = 100 

min 

times  the  previous  values.  Thus,  elimination  of  K /K  • 

^ ’ max  min 

in  highly  desirable  and  the  SqAL  (Fig.  4a)  is  a 

suitable  structure  for  this  purpose. 


III.  SOAL  ADAPTIVE  THEORY 


In  order  to  eliminate  K /K  . from  Yt  in  (9), 

max  min  '1  ’ 

it  is  necessary  to  eliminate  K from  (8)  and  K . 

from  in  the  discussion  leading  to  (9).  Both  are  due 

to  K/A  being  invariant  in  (7b),  which  in  turn  gives  the 

system  zero  sensitivity  to  uncertainty  in  K . Hence, 

K elimination  with  retention  of  zero  K-sensitivity 

max  min  ^ 

is  possible,  if  instead  of  K/A  , BK  in  (7b)  is  made  ■ 
invariant.  But  whereas  A automatically  tracks  K such 
that  in  steady-state  K/A  is  invariant,  B can  be  made 
to  do  so  only  by  means  of  a secondary  loop.  The  function 
of  the  latter  is  to  track  A and  adjust  B , to  force  A 
to  be  constant  in  the  steady-state.  Since  (7a)  fixes  w 

o 

and  (7b)  compels  BK/A  to  be  fixed  (inasmucii  as 

fixed^  fixed  A compels  BK  to  be  invariant, 
This  rcasining  leads  directly  to  the  self-oscillating 
adaptive  loop  (SqAL)  structure  of  Fig.  4a. 
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SOAL  Structure.  .In  Fig.  4a,  a peak  detector  A is  used  to 

measure  the  difference  between  the  maximum  and  minimum 

values  of  x =A  sin  w^t  + x^(t).  This  gives  a fast  and 

accurate  measurement  of  2A,  because  (recall  2a,  b)  the 

system  is  designed  so  that  over  the  range  of  command  and 

disturbance  inputs,  jx^Ct) [ is  both  small  and  slow 

relative  to  A and  . (See  later  for  comparison  vjith 

other  methods).  The  measured  A^  is  compared  with  the 

chosen  reference  value  A and  the  difference  A (t)  , 

r e 

processed  by  )/'(s)  , changes  B so  that  in  the  steady- 
state,  A^  is  very  small  or  zero.  The  dynamics  of  the 
secondary  loop  are  next  considered. 


Dynamics  of  Secondary  Loop 


The  technique  [293  of  the  'Fundamental  Feedback  Equation' 
is  used  to  derive  the  dynamics.  In  a linear  time  invariant 
system,  let  be  any  variable  chosen  cis  an  independent 
system  input,  and  O'  the  output  variable  of  interest.  Let 
(s)  be  the  transfer  function  of  any  system  element  such 
that  /^(s)  = JSCs)  ( s)  , where  J>  and  /e  are  any  two 
. system  variable  so  related  to  ^ . If  the  element  is 

thus  uniquely  identified  by  the  variables  and  ^ , 

then  the  relation  between  ^ and  , has  the  form 


= -A/rr  = t_.(s)  + — . 


(11) 


1 


The  important  feature  of  this  'fundamental  feedback 
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equation'  for  reference  ^ (s)  , is  that  t^^  , t^^  , t^^  , 

t (s)  are  completely  independent  of  /^(s)  and  they  are 
calculated  quite  easily,  as  will  be  seen.  Eq.  (11)  is  there- 
fore very  useful  if  one  wishes  to  focus  attention  on  the 
relation  between  and  It  is  also  very  useful  in 

analysis  [29]  enabling  a complex  problem  to  be  broken  into 
four  simpler  ones-so  used  here. 

Let  the  independent  input  ^ be  a change  in  A and 
the  output  ^ be  the  change  in  B . Let  ^(s)  = (s)  , 

so  A s)  and  (s)  are  as  shown  in  Fg.  4^1.  The  four 
transmissions  t^j(s)  are  obtained  as  follows  [29].  Cut 
the  system  at  points  a,  a'  with  Fig.  4b  replacing  the  top 
(secondary  loop)  part  of  Fig.  4a.  The  functions  > ^ci  the 

transmissions  from  = AA  to  the  variables  (9-  =AB  , 

^ = a . The  functions  t^^  , t^^  are  the  transmissions 
from  ^ = a'  regarded  as  an  independent  inpu-^  to  the 
variables  =AB  , ^ = a . One  must  calculate  four 
transfer  functions  instead  of  only  one  - but  the  cuts  can 
make  the  latter  four  much  easier  to  find  - providing  a wise 
choice  of  has  been  made.  The  system  input  r(t)  in 

Fig.  4b  need  not  be  zero  so  long  as  it  may  be  assumed  it 
does  not  affect  the  limit  cycle. 


Accordingly,  in  Fig.  4b 

a 


£ 


^ , = 0 ; t . = -r-l  - -A,  ^(s)  being 

/I  Z^-0 

the  transfer  function  of  a peak  - to  - peak  detector,  which 
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requires  a half  - period,  giving 


'Cl 


-e 


-.5sT 


T = 2Tr/a) 


os 


The  function  t 


sc 


't?=0 

consists  of  two  parts. 


= 1 


= 0 


Obviously  one  part  is  -A  (s)  . The  other  is  the  effect 
^ me 


of  AB  on  the  limit  cycle  amplitude  at  x , precisely  as 
in  the  SOAS.  The  steady-state  value  is  obtained  from  (7b) 
giving  ^ ^ I ' ‘ dynamics  of  the 

f 

change  AA  due  to  AB  was  found  experimentally  to  be 

Lf/P 

well  approximated’  by  with  p between  2 and  3 

for  different  conditions.  The  worst  case  P = 3 is  used 
here,  giving  the  ’secondary  loop  transmission' 


L^(.) 


-*tcs  ' e 


-.5sT 


L^/3 

l+L^/3 


(12) 


The  stability  and  speed  of  this  B adjusting  loop  is 

primarily  determined  by  L (s)  , whose  maximum  'crossover 

I I ^ 

frequency'  (where  |Lg(j<^j2)|=  1 ) is  < w^/3 

— 5sT 

approximately,  because  e ' contributes  -60°  phase 

“o 

at  — 2 and  allowing  at  least  30°  phase  margin,  thefe  is 

left  -90°  needed  for  |dLg( ja))/dm  [;«  -6db  per  octave 
at  oj^  . It  is  interesting  to  see  the  reappearance  of 
the  K uncertainty  in  , after  it  was  banished  from 

the  primary  loop  and  quasiline^irity  conditions.  The  K 
uncertainty  must  be  considered  in  the  shaping  of  L 


Since  w 


w(o/3  at  K-K  , w • at  K = K. 
c max  o max  cmin  min 
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can  be  quite  small  with  resulting  much  slower  adjustment  of 

B to  K changes  near  K = K . . The  function  <I)(s)  is 

^ min  ^ 

available  for  shaping  of  Lg^s)  . 


From  the  above  it  is  easy  to  find  the  effect  of  changes 
in  K on  B . The  independent  input  ^ is  now  a change  in 
K . Using  the  same  A (s)  , the  result  is 


AB  B ^s 

AK  ■ " K 1+Lg 


L of  Eq.  (12). 
s 


SOAL  Design  Equations 

These  follow  easily  from  those  of  the  SOAS.  In  (8)  BK 

is  constant,  so  K is  replaced  by  K . Constraint  (9) 

' max  ^ 

becomes 


[X  (jw  )|  _ 

! -A---  > h 


m e o 


(13) 


because  A is  constant.  For  the  same  reason  (10)  becomes 


A~  ^ 5^’  » V 0.  .<  . 


(14) 


The  above  is  illustrated  by  a numerical  example. 


SOAL  Design  Example. 


The  plant  and  specifications  are  the  same  as  in  the 
SOAS  example,  except  that  m there  of  1 is  decreased  here 
to  .01.  In  this  way, 


K 


max 


— ^ |Z  (jm  )|  (SoAL)  = — 
'e-’o' 
new  old  min 


K 


|Z  (jio  ) 1 (SOAS) 


giving  again  =4.5  [If  m = .01  is  used  in  the 

SOAS,  the  resulting  to^  = 14  rps , and  it  is  possible 

that  over  part  of  this  range"  |Z^^(ja))|  > |Z^g(ja))|  , so 

even  larger  would  be  necessary].  The  bounds  on 

L^Cju)  to  obtain  the  specified  disturbance  attenuation 

are  the  same  as  in  the  SOAS  so  precisely  the  same  L^(s) 

Go(s)  , BG^Cs)  , G-,Cs)  result.  In  the  SOAS,  A . = . 5ti 

■32’!  » min 

here  fixed.  If  B = IOtt  is  used  at  K • = 1 , then 

min  ' 

B = .16ir  at  K . = 100  . Only  ip(s)  remain  to  be  chosen 

mdx 

Choice  of  tp(s)  = 0.4/s  gives  very  simple  secondary  loop 
dynamics,  because  at  K = 100  , A = .5  , B = lOir  , 

Lg  » I ^ (ignoring  (L^/ 3)/[l  )]  which  is 

close  to  1 over  [0,.4]  , .4rps  being  the  crossover 

frequency  m of  L (ju)  . The  time  del^iy  e~*^^  gives 
c s 

16®  phase  lag  at  . If  this  too  is  ignored. 
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Design  Verification.  The  system  was  simulated  on  the 
CDC  6400  at  University  of  Colorado,  using  MIMIC.  Results 
are  shown  in  Figs.  5a-c  for  abrupt  gain  changes  in  K . In 
(a),  the  time  required  for  B to  reach  95%  of  its  final 
value  is  9.8  seconds.  If  K = 100  is  used  in  the  previous 
simple  first  order  approximation  of  L then  three  time 
constants  plus  a half-period  predict  8.2  seconds,  whereas 
use  of  K = 50  predicts  15.7  sec.  It  is  seen  that  the 
speed  of  adjustment  significantly  decreases , as  predicted, 
at  the  smaller  K values.  The  invariance  of  the  system 
step  response  for  various  K is  seen  in  Figs.  6a-^-,  also 
that  (2a)  is  very  well  satisfied. 

/ 

Comparison  with  other  Two-Loop  Self-Oscillating  Adaptive 
Systems . 

One  significant  distinction  of  the  SOAL  is  its  retention 
of  the  nonlinear  element  N of  the  SOAS,  with  use  of  the 
second  loop  to  modify  its  parameter.  In  some  other  two-loop 
systems  [13,  18,  19,  2lj , N is  eliminated  and  a gain 
changer  is  used  in  its  place.  Another  crucial  difference  is 
the  point  at  which  A of  Asinu^t  is  measured.  In  the 
SOAL,  it  is  measured  at  the  point  where  it  has  the  greatest 
immnunity  to  forced  signal  components,  because  the  design 
guarantees  that  at  x in  Figs.  1,  4 the  forced  components 
are  slow  and  small  realtive  to  and  A . Peak  to  peak 

detection  becomes  possible  which  further  increases  the 
immunity.  In  the  others  [18-221,  A is  measured  closer 
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to  the  plant  input  or  output,  where  the  forced  components 
are  much  stronger.  A bandpass  filter  becomes  necessary  vjith 
its  additional  significant  time  log  over  and  above  that 
needed  to  measure  A . Disturbance  components  can  conceivably 
be  misinterpreted  as  plant  gain  changes  v;ith  possibly  disastrous 
results.  If  the  pass-band  of  the  filter  is  narrowed  to 
decrease  the  latter  possibility,  the  greater  is  the  danger 
if  wanders  out  of  the  pass-band  and  the  greater  is  the 

time  lag  of  the  filter. 


General  Plant  Uncertainty. 


In  the  above,  only  uncertainty  in  the  plant  gain  factor 
K of  P = KP^  was  considered,  in  order  to  emphasize  the 
fundamental  SOAS  and  SOAL  properties.  In  the  general 
uncertainty  case  has  adaptive  obligations,  which  it 

does  not  have  for  K variations  only  and  no  disturbances. 
These  obligations  in  a quantitative  design  can  be  put  in  the 
form  of  bounds  on  L^Cjoj)  , just  like  those  in  Fig.  3 for 
the  disturbance  attenuation  obligations  of  .' 


Reference  1 25]  has  treated  in  detail  the  general 
uncertainty  problem  for  the  SOAS  and  it  was  found  that 
some  iteration  may  be  theoretically  necessary,  but  rarely 
in  practice  because  is  relatively  so  large  that 

|L^(jw)l  is  significant  over  a relatively  large  w range. 
The  treatment  for  the  SOAL  is  precisely  the  same  and  there- 
fore need  not  be  repeated  here.  However,  in  the  design 
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execution,  the  iteration  is  much  more  likely  to  he  necessary, 

especially  in  the  case  of  large  uncertainty  in  plant  dyruimics. 

The  reason  is  that  u is  much  smaller  in  the  SOAl,  than  in 

o 

the  SOAS  5 so  one  does  not  have  so  much  extra  bandwidth. 

The  bounds  on  due  to  the  uncertainty  in  1^ 

be  such  that  needed  (at  vjhich  = -0.5)  for  this 

purpose,  is  larger  than  o)^  min  obtained  from  constraints 
(13,  14).  It  is  obviously  easy  to  manufacture  such  problem 
specifications.  In  any  case,  the  procedure  is  identical  to 
that  presented  [25]  for  the  SOAS. 

IV.  NONLINEAR  SECONDARY  LOOP  SOANL 

The  logical  next  step  for  improving  the  addiptive  system 

is  to  try  to  eliminate  the  sensitivity  of  the  SOAL  secondary 

loop  Lg  to  the  plant  gain  K . This  can  be  done , in  the 

steady-state,  by  means  of  a multiplier  in  the  secondary  loop 

whose  multiplication  factor  is  inversely  proportioniil  to  K . 

Since  K is  not  directly  available,  B which  is  inversely 

related  to  K is  used,  as  shown  in  Fig.  7,  which  shows  only 

the  secondary  loop,  as  the  primary  loop  is  identical  to  that 

in  Figs.  1,  4a.  The  ’fundamental  feedback  equation'  is 

again  used  to  find  the  dynamics  of  the  new  L^  . As  the  loop 

is  nonlinear,  such  a linear  time  invariant  description, 

valid  only  for  small  perturbations  about  an  operating  point, 

is  a function  of  the  operating  point.  The  easiest  one  is 

when  the  secondary  loop  is  initially  quiescent  i.e. 

A = A so  there  is  no  multiplier  output. 
r in  t'  t 


1 


I 
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Let  ^ and  ^ be  as  shown  in  Fig.  7.  Let  ^ be 
AK  , a change  in  the  plant  gain,  and  - Ab  . Make 


cuts  at  a , a*  . Then  t . =0  due  to  the  cuts  and  t =1 

oi  os 


To  find  t . = XT?  > note  that  AA  (recall  similar 
Cl  AK  m 


calculation  in  the  SOAL) 


4B(AK) 


L.-/3 

' 2 h 1 o ' me  l + L^/3 


This  is  multiplied  by  ^ result  proportional 


to  B . To  find  t^^  , let  a'  be  an  independent  input 


and  find  a/a'  . B becomes  B = B + ^ 

new 


giving 


AA,„  = | G„P,  G,  ( jw^)  1 K A 

m TT  ^ n I r\  n 


L^/3 


which  is  multiplied  by  -B  /B  . . For  incremental  a' 

new  min 


the  result  is 


B . a' 
min 


, and  L = -ip  (s)  t 

s . sc 


is  the  SOANL  secondary  loop  transmission.  In  the  steady 


state  the  above  calculated  AB/AK(s)  gives  AB/AK  = -B/K  , 


just  as  in  the  SOAL  . The  important  difference  is  that  L 


contains  BK  (cf  K in  (12)  for  the  SOAL) , which  is 


invariant  in  the  steady-state.  Hence,  the  SCAIIL  should 


deliver  for  all  K G [K  . , K ],  the  foist  adjustment  to 

min  ’ max  ’ 


K variations,  which  the  SOAL  can  do  only  near  K 
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SOANL  Dos  ip, n Example 

The  synthesis  procedure,  specifications  etc  are  the 
same  as  for  the  SOAL.  Hence,  the  results  for  , G2 , B 

etc  are  precisely  the  same.  Simulation  results  are  shown 
in  Fig.  8.  The  analysis  made  was  for  small  K changes, 
and  would  be  expected  to  be  valid  for  large  changes  if  these 
were  somewhat  slower  than  the  secondary  loop  dynamics. 
Nevertheless,  the  same  basisc  results  hold  even  for  very 
large  and  abrupt  changes  in  the  plant  gain  K , as  seen  in 
Fig.  8. 

V.  CONCLUSIONS 

The  development  of  a quantitative  design  theory  for  the 

single  loop  self-oscillating  adaptive  system  (SOAS)  has  led 

in  a natural,  logical  manner  to  the  two-loop  system  (SOAL)  v;ith 

important  superior  features.  The  emphasis  on  quantitative 

design  reveals  also  the  superiority  of  the  SOAL  over  other 

such  two- loop  structures.  It  is  interesting  that  the  gain 

was 

factor  uncertainty  ^max'^^min  ^Bichyyinherently  banished  as 

an  adaptive  problem  by  the  oscillation,  returns  in  the 

SOAS.  in  the  quasilinearity  constraints.  It  f)lays  there  an 

important  role  in  forcing  an  unnecessarily  large  . 

It  is  banished  in  the  SOAL  permitting  a mucli  smallei"'  to  . , 

but  reappears  in  the  secondary  ]oop  dynamics.  This  results 

in  the  speed  of  adjustment  to  K variations  being  a function 

of  K and  much  slower  at  K • tlian  at  K 

ii\  1 n 


max 


It  is 


1 
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banisTicd  from  the  SOAL  by  moans  of  nonlinear  compensation 

in  the  secondary  loop,  giving  the  SOAh'L..  It 

undoubtedly  reappears  in  a subtle  manner  in  a nonlinear 

analysis  of  the  SOAIJL  > but  if  so,  it  does  not  appear  to 

be  a critical  one.  Exactly  the  same  tv;o-loop  extensions 

can  be  made  to  the  externally  excited  adaptive  (EEAS) 

'•  EEAL  and 

system  [26],  giving  multiple  loop^EEANL,  V'/hich  therefore  do  not 

require  separate  treatment. 
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Design  constraints  on  Bode  plot 
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Fig.  5.  SOAL  Desi^ 
Example, 

Effect  of  abrupt  K 
chan^ges  on  B(t)  and 
X It)  = A(t)  sino)  t. 

(a)  K initial  = 100 
final  = 50 


(b)  K,  initial  = 50 


Invariance  of  step 
response 

(a)  at  K = 100 


